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Abstract 

We show that the gauge and metric field contribution to the axial anomaly of a four- 
dimensional massless Lifshitz fermion theory with anisotropy scaling exponent z is 
identical to the relativistic case, hereby extending the results found in arXiv: 1103. 5693,, 
to arbitrary values of z. This is in accordance with the fact that the axial anomaly is 
an infra-red phenomenon in disguise. We also provide some new models that realize 
baryon and lepton number violation in non-relativistic theories of gravity. In all cases, 
the volume of space exhibits a lower bound that is fixed by the gravitational coupling 
parameters. 



1 Introduction 



It has been known for a long time that a massless Dirac fermion in four space-time 
dimensions exhibits a quantum anomaly in the axial current conservation law when 
coupled to external gauge and/ or metric fields. The anomalous term in the divergence 
of the axial current = Yjl-^j^W is proportional to the characteristic classes Tr(f A f ) 
and Tr(J^ A -R) written in terms of the curvature two-forms F and R"jj of the gauge 
and metric fields, respectively, with corresponding coefficients l/47r^ and l/967r^ (in 
appropriate units that will be used throughout this paper) IH |2l. The coefficient of 
anomaly is half of that when computed for a Weyl fermion. This result describes the 
local form of the axial anomaly, whereas its integral over space-time has been used 
to study possible violations of chiral symmetry via the Atiyah-Singer index theorem 
(for an overview of the subject from a physicist's perspective see, for instance, [3] and 
references therein). The axial charge Q5 is not conserved when the four-dimensional 
Dirac operator in a given background exhibits an asymmetry in the number of positive 
and negative chirality zero modes, leading to violations of lepton and baryon number 
in the theory BU. Generalizations to higher dimensions have also been worked out 
systematically [5] and the results found many important applications in physics. 

In this paper we examine the occurrence of axial anomalies in the non-relativistic 
fermion field theories of Lifshitz type. For this we consider higher order generaliza- 
tions of the Dirac operator of the form Q = ij^do + iYdi{—dj(d^ + M^)* acting on 
spinors Y, where M is an arbitrary mass scale and a is a free parameter. The result- 
ing theory of Lifshitz fermions, which is described by the Lagrangian density YQY 
exhibits anisotropic scaling in space and time by assigning dimensions [L] = — 1, 
[T] = —z and [Y] = z/2 with parameter z = 2a + 1. This model has an axial current 
which is conserved at the classical level. Its time component = Y7'^75Y is the same 
as in the relativistic case (the same also applies to the axial charge Q5), but the spa- 
tial components are more complicated as they involve derivative terms that depend 
upon a; their precise form will not be needed for the purposes of the present work. 
Then, a natural question is whether there are anomalies in the axial current conser- 
vation law that arise upon quantization of the fermions in the background of gauge 
and/ or metric fields, and, if so, what will be their form compared to the relativistic 
case, a. = 0. The result turns out to be independent of the parameters M and a, and, 
in this sense, the structure of the axial anomalies appears to be universal. 

The present work generalizes earlier results on the subject that were obtained for 
Lifshitz fermions with anisotropy scaling parameter z = 3 coupled to background 
gauge fields [6] and gravity [7J in four space-time dimensions. Here, z = 2a; + 1 is left 
arbitrary and it can even assume fractional values. For practical reasons we restrict 
attention to four space-time dimensions, although we expect the same universality to 
govern the form of the axial anomalies in higher dimensional Lifshitz fermion models 
as in relativistic theories |]5|]. An intuitive explanation for all this is provided by the 
fact that the axial anomaly is an infra-red phenomenon in disguise, and, hence, any 
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higher order corrections to the Dirac operator that become relevant in the ultra-violet 
regime ought to leave its form unaltered (the way it works in relativistic theories is 
nicely explained in ||5]], but see also ||8l HI [101 ^or earlier work on the subject). The 
introduction of the mass scale M into the definition of Q serves precisely the purpose 
of showing by explicit computation that the axial anomaly is indeed independent of it 
and the order of the higher derivative corrections to the fermion operator. The result 
should be contrasted to the form of the Weyl anomaly in Lifshitz theories that exhibit 
anisotropic local scale invariance. In that case one does not expect the result to be the 
same as for relativistic fields because the Weyl anomaly is an ultra-violet phenomenon 
that is sensitive to higher order corrections at short distances. Some partial results in 
this direction can be found in |11J for scalar field models, whereas generalizations to 
other field theories of Lifshitz type have not been addressed so far in the literature. 

Thus, the computation of axial anomalies in non-relativistic massless fermion mod- 
els is a much simpler problem compared to other type of anomalies and at the same 
time it has many interesting applications related to chiral symmetry breaking at a 
Lifshitz point. A closely related problem is the parity anomaly of three-dimensional 
fermion theories coupled to external gauge and/ or metric fields. The methods we use 
in this paper suffice to demonstrate the universality of the parity violating piece of 
the effective action induced by massless fermions regardless the order of the three- 
dimensional fermion operator. This result is a manifestation of the close relation that 
exists between different type of anomalies in odd and even number of space-time 
dimensions generalizing the framework that was established long time ago for rela- 
tivistic theories (see, for instance, ||T2l| and references therein). We will say more about 
this later focusing, in particular, to the gravitational case and the universal character 
of the induced three-dimensional Chem-Simons action. 

In section 2, we employ the path integral method to extract the local form of the 
axial anomaly of Lifshitz fermion models. The computation is performed rather eas- 
ily for gauge field backgrounds and the result is identical to the relativistic case, as 
expected on general grounds. Confirming the universality of the axial anomaly for 
metric field backgrounds is much more involved computationally by path integral 
methods. For this reason we employ an indirect method based on the integrated form 
of the anomaly and the associated index theorem for the Dirac-Lif shitz operator. Then, 
using this framework, the coefficient of the axial anomaly follows from the computa- 
tion of the ;y-invariant of the associated three-dimensional Lifshitz operator on certain 
geometries. In section 3, we briefly review the notion of ^/-invariant and its relation 
to anomalies in relativistic fermion theories in three and four dimensions. In section 
4, we compute the //-invariant of higher order fermion operators on homogeneous 
geometries and show that the result is independent of the parameters M and a. This 
suffices to prove our claim about the gravitational contribution to the axial anomaly of 
Lifshitz fermion. In section 5, we present some applications to gravitational theories of 
Lifshitz type and examine the possibility of chiral symmetry breaking by instantons, 
thus generalizing the results found in [7J. Finally, section 6 contains our conclusions. 
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2 Path integral derivation of axial anomaly 



Coupling the Lifshitz fermion theory to external fields amounts to replacing dpi by Dpi, 
as usual. Then, the interacting Dirac-Lifshitz operator takes the form 

Q = i^Vp, = h^Do + izV[D,(-DfcD'^ + M^r + {-D^D^ + M^fDi] (2.1) 

up to a factor ordering ambiguity that turns out to be irrelevant for the computa- 
tion of the anomaly. To avoid unnecessary complications or ambiguities in the def- 
inition of the operator (|2.1|l . and in view of the applications that will be discussed 
later, coupling to geometry is taken with respect to metrics — {dx^Y + Siji.^'^' x)dx^dxK 
Such backgrounds arise naturally in the canonical decomposition of the metric field 
in space-times of the form I x E3 choosing the lapse and shift functions as N = 1 and 
Ni = 0. 

The anomalous divergence of the axial current can be found by applying Noether 's 
procedure to the fermionic path integral of the interacting theory, which is most con- 
veniently described in the Euclidean domain after Wick rotation of time x^ — > it. The 
measure (PY) (PY) is not invariant under chiral rotations = 26:75 Y, giving rise 
to a non-trivial contribution to the divergence of the axial current (see, for instance. 



| |T3| and references therein to the original papers). Using the eigen-states (pn{t,x) of 



the interacting fermion operator i^y^^Dpi, the result takes the intermediate form 

V,/f (f,x) = 2 lim Zcplit,x)j5 e-(''^"^/')'/^'>„(i,x) (2.2) 



after introducing a cut-off A to regulate the infinite sum that otherwise is ill-defined. 
Note that A is raised to the power 2z to match the scaling properties of the square of 
the Dirac-Lifshitz operator for general values of z. 

The gauge field contribution to the axial anomaly can be computed relatively easy 
for all values of the anisotropy scaling parameter z. First, using the algebra of the Dirac 
gamma-matrices, it is convenient to rewrite the square of the Dirac-Lifshitz operator 
appearing in the regulator as 

ih^Vpi)^ = -V,V^ - 7l [D,, Vy] , (2.3) 

where = dpi — iApi provides the minimal coupling of the theory to an external 
gauge field A^, which can be Abelian or non-Abelian. Next, using the plane wave 
basis of solutions of the free Dirac-Lifshitz operator, exp(zfc^x^), we can extract the 
explicit gauge field dependence of the anomalous term that appears on the right-hand 
side of equation (|2.2|l . Since the action of the interacting Dirac-Lifshitz operator on 
plane waves amounts to replacing by Dpi + ikpi everywhere, it follows that the 
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anomalous divergence of the axial current in a gauge field background takes the form, 
V,Jl = 2 km^A^« Tr/ ^ .-"l-"^ exp{ - (f„, + 2.F„,^.*') 

- ^ h', 7"] (f/t + i'xFiM') +■■■}. (2.4) 

In writing this equation we have also rescaled Rq to A^Rq and kj to Afc, conforming 
with the anisotropic scaling of the components of the momentum vector = {ko, ki) 
that is implied by the scaling of the space-time coordinates = [xq, xi) and intro- 
duced the unit three-momentum vector with components ki, setting ki = kki, where 
k^ = kikK Here, we have depicted the most relevant terms that originate from the 
exponential of {i^^V^iY after replacing D^i by D^i + ik^. The factor exp(— /cq — k^^) 
originates from V^V^ in the decomposition (|2.3|) , which is taken in the Euclidean do- 
main, whereas the two terms shown explicitly in curly brackets originate from the 
electric and magnetic components [7*^, 7'][Po/ T^i] and [7^, J^j^Dj, "Dj^], respectively, 
appearing on the right-hand side of (|2.3|) . The terms that are omitted carry explicit 
dependence on the mass scale M, but they do not contribute to the final result. Such 
terms yield contributions that either vanish as A ^ 00 or their trace is zero because 
they do not contain sufficient number of gamma-matrices multiplied with 75. 

The calculation proceed by expanding the exponential in curly brackets and noting 
that the only relevant terms arise to second order. We also perform the Gaussian inte- 
gration over kQ, picking up a factor of i/tt, and use the identity Tr(75 [7*^, 7'] [ji, j^] ) = 
—16e^^^^. Then, after taking the limit A — 00, we obtain 

^?^^5 = / d^ke-^'^k'"~^^r (Fq^F^-, + 2aFoiFjM' + ^ocForFjM) , (2.5) 

setting z = 2a: + 1. The calculation is completed by integrating out the three-momenta. 
Using polar coordinates in momentum space with k = (sin^sin^, sin0cos(/), cos9) and 
substituting x = k^, it follows easily that 

J e-^'^k^'-^ = 471 JJ' dk k^^-h-^"' = ^7r3/2 (2.6) 

and 

j d^k e-^'^k^--^ kkj = ^ dk k^'-'e-''' = . (2.7) 

The final result is independent of the scaling parameter z, 

V./f = 4^^°'^'Tr(Fo,F^O = 4^Tr(F A F) , (2.8) 
and coincides with the anomaly of the axial current conservation law of a Dirac f ermion. 
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Similarly, we may compute the gravitational contribution to the axial anomaly un- 
der the minimal substitution of 9^ by D^j = 9^ + (1/8) [7^, 'yi,]co^^^ written in terms 
of the spin connection OJ. The regulated sum shown in (|2.2|) is much more difficult to 
evaluate in this case as one has to introduce point splitting and use the curved space 
analogue of plane waves given in terms of the geodesic interval in space-time. Keep- 
ing track of the relevant terms is also quite tricky, in general, apart from the z = 1 
case that corresponds to Dirac fermions. For higher values of z one encounters mul- 
tiple derivatives of the geodesic interval which in the coincidence limit correspond to 
Synge-DeWitt tensors of rank bigger than 2. There are many sources of such terms 
that contribute to the final result and one has to sum them up carefully For z = 3 
the intermediate calculations are already quite tedious (see [Zl for details) and one ex- 
pects that this method will be practically impossible to carry out to the end for general 
values of z. A closely related (but practically simpler) calculational method could be 
developed by applying path integral methods to suitably chosen models of supersym- 
metric quantum mechanics with Hamiltonian equal to the square of the Dirac-Lif shitz 
fermion operator. This method was successfully employed in the past to compute the 
form of anomalies in relativistic theories ||5]], particularly in higher dimensions where 
the conventional approach becomes intractable, but it is not currently known how to 
generalize it to encompass non-relativistic fermion systems. Thus, we should resort to 
other means for the efficient calculation of the gravitational contribution to the axial 
anomaly for general values of the anisotropy scaling parameter z. 

The axial anomaly is necessarily a total derivative term, as it obstructs the diver- 
gence of a current. It is also gauge invariant, meaning that the anomalous term is 
solely expressed via the field strength F and/ or R. In four dimensions there is only one 
such term given by the topological density Tr(F A F) when fermions couple to gauge 
fields. Thus, the whole purpose of the path integral computation is to determine the 
coefficient of this term, which, as we saw before, turns out to be universal. By the same 
token, when fermions couple to metric fields, there are two possible terms given by the 
topological densities Tr(R ^R) = R% A R^a and Tr{R A *R) = {-l/2)eabcdR% A R'd 
in four dimensions. The second term is easily ruled out by the trace identities of 
Dirac gamma-matrices; close inspection of the anomalous term (|2.2[) reveals that the 
trace cannot yield more that one epsilon symbol, and, hence only Tr{R A R) and not 
Tr(R A*R) can possibly contribute to the axial anomaly. Then, the only remaining task 
is to evaluate its coefficient and verify that in all cases the result is the same, namely 

V./f = - ^2^''''R''\rRa, = A R) . (2.9) 

The method that we will follow relies on the integrated form of the anomaly and its 
relation to the index theorem of the Dirac-Lifshitz operator. An important ingredient 
in the computation is the ^/-invariant of the fermion operator restricted on the three- 
dimensional leaves of space-time foliation and its intimate relation to the gravitational 
Chern-Simons action via the anomaly, which turns out to be universal. 
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3 The //-invariant and axial anomalies 



In this section we make a small detour to summarize some basic facts about the index 
theorem of the Dirac operator and the associated notion of //-invariant, which play 
important role in relativistic quantum field theory. Later we will extend the results to 
the more general class of Dirac-Lifshitz fermion operators and use them to determine 
the coefficient of the axial anomaly by spectral methods. 

Recall that the 7/-invariant of a self-adjoint matrix A with real eigen-values Aj^ is the 
signature of A, which is defined to be 

y]A = sign(A) = E '(sign A^) = E ^ " E ^ ' (3-1) 

assuming that there are no zero eigen-values in the spectrum, which otherwise are ex- 
cluded. As such, r]j{ is invariant under rescaling of the matrix by any positive number 
and it counts the spectral asymmetry between positive and negative eigen-values of 
A. If A depends upon a parameter S, the spectrum will flow as S varies continuously 
and level crossing may occur in the system. The //-invariant of A(S) changes by ±2 
units every time a negative (respectively positive) eigen-value at crosses to positive 
(respectively negative) values at Si as S varies from Sq to ^i. 

Computing the //-invariant is a notoriously difficult task when the matrix A is in- 
finite dimensional, since the spectrum of A cannot be determined in general and the 
definition of //^ becomes ill-defined as it stands. The standard prescription in those 
cases is to consider first the //-function of A obtained by zeta-function regularization, 
as 

7a(s) = E '(sign A|c)|Afc^^ (3.2) 

h 

and then define the corresponding //-invariant as 

Va = Va{0). (3.3) 

Typically, A is a self-adjoint operator on a compact manifold so that its spectrum is 
discrete and //^(s) can be analytically extended to a meromorphic function on the 
entire complex s-plane which is regular at s = 0. Yet, direct computation of the //- 
invariant by spectral methods can only be done in a few cases. Such examples will be 
given later. 

The Atiyah-Patodi-Singer index theorem is the main mathematical framework that 
necessitated the introduction of the //-invariant of various operators fill (but see also 
the textbooks IITSi 1161 for more details). Restricting attention to the Dirac operator 
/7^'D^ on a four-dimensional Riemannian manifold M4, the difference of its positive 
and negative chirality zero modes, known as the index of i^^D^ (we will use D in 
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short), is given by 



Ind(D) = / Tr(R A R) + / Tr{e A R) - U^{dM^) . (3.4) 

The first term provides the bulk contribution to the index, which is obtained by in- 
tegration over M4 of (a half times) the local form of the axial anomaly of a Dirac 
fermion in the background of a metric field. It involves the (appropriately normal- 
ized) characteristic class Tr{R A R) that is written in terms of the curvature two-form 
R = dco + CO A CO associated to the spin-connection co of the metric g on M4. The other 
two terms arise when M4 has a boundary 9M4. The local boundary term is the integral 
over 9M4 of the (appropriately normalized) secondary characteristic class Tr{6 A R), 
lUZl, which is written in terms of the second fundamental form 6 = co — co^ and ac- 
counts for the possible deviation of the metric g from a cross-product form gQ at the 
boundary. Finally, the last term is a non-local boundary contribution to the index pro- 
vided by the ;y-invariant of the tangential part of the Dirac operator restricted to 9M4. 
Here, it is implicitly assumed that the four-component spinors Y have a specific fall- 
off rate close to dM^ (known as APS boundary conditions) and so what we are really 
computing is the L^-index of normalizable zero modes of iy^D^. 

An immediate consequence of the Atiyah-Patodi-Singer index theorem is that the 
//-invariant of the three-dimensional Dirac operator on a compact manifold E3 without 
boundary, which is endowed with a Riemannian metric, is related to the gravitational 



Chern-Simons action |T7| 



Wcs(2^3) = J 1r (^co Adco +'^co Aco Aco^ (3.5) 

as follows, 

lf]D{^3) + T^^cs{^3) = constant . (3.6) 

This is easily seen by considering a four-dimensional spin manifold M4 bounded by E3 
and choosing a metric on it that has cross-product form near the boundary dM^ = L3. 
Thus, although both the //-invariant and the Chem-Simons action depend explicitly 
on the conformal class of the metric on E3, their (appropriately weighted) sum is in 
fact independent of it. The coefficient of Wcs is uniquely fixed by the anomaly. We 
also note for completeness that both i/d and Wcs ^re odd under parity, thus flipping 
sign under orientation reversing diffeomorphisms on the manifold E3. 

Equation (|3.6|) provides an indirect way to compute the t] invariant of the three- 



dimensional Dirac operator, up to a constant, based on axial anomalies in four space- 
time dimensions. Conversely, one may use this formula to extract the coefficient of 
the anomalous term Tr{R A R) in the divergence of the axial current using spectral 
methods, provided that the //-invariant can be computed explicitly for a certain class 
of three-metrics. This is the strategy that will also be applied later to Lifshitz models. 



7 



Alternatively, on a four-dimensional manifold of the form M4 ~ 7 x E3, where I is 
an interval (it can also be J ~ R in the physical applications) endowed with the metric 

ds^ = dt^ + gij{t, x)dx^dx^ , (3.7) 

the Atiyah-Patodi-Singer index theorem for the Dirac operator implies the following 
useful relation 

^A7d(5:3) + ^^AWcs(i:3) = integer , (3.8) 

where A denotes the difference of any given quantity at the two end-points of the 
Euclidean time interval I. One may view (|3.6|) as special case of the more general 
relation (|3.8)) by just shrinking one of the two ends of the cylinder J x Z3 to a point. 
For example, for E3 ~ S^, let us assume without loss of generality that one of the two 
boundaries is endowed with the constant curvature (round) metric. At that boundary 
the //-invariant vanishes, as there can be no asymmetry in the spectrum of the Dirac 
operator, whereas W^s assumes the value 167r^. Shrinking that boundary to a point 
by conformal rescaltng does not affect the values of ?/d and W^g and (|3.8)) reduces to 
(|3.6|) with constant term equal to 167r^/1927r^ = 1/12 modulo integers. 

The integer appearing on the right-hand side of equation (|3.8|) is the index of the 
Dirac operator on J x E3. If the space-time metric (|3.7|) has a cross-product form every- 
where, i.e., gij is independent of t, the three-dimensional metric at the two ends of the 
cylinder (boundaries of space-time) will be the same, and, hence, A//d = = AW^s 
giving zero index. Otherwise, if the metric is warped by depending explicitly on f , the 
//-invariant and the Chem-Simons action will be different at the two end-points in gen- 
eral. Furthermore, it can be shown in this case that the index of the four-dimensional 
Dirac operator is provided by the spectral flow, 

Ind(D) = AS(E3) , (3.9) 

where AS(E3) is the net number of level crossings that may occur in the spectrum of 
the three-dimensional Dirac operator on E3 as the metric gij{t,x) deforms from one 
end-point of the interval I to the other end-point (recall that the //-invariant jumps by 
±2 units when an eigen-value crosses from negative to positive values or conversely, 
and, therefore, each level crossing contributes ±1 units to the index). This formula 
and its generalization to Lif shitz models is very useful for the applications that will be 
discussed later. 

There is another occurrence of the //-invariant of the three-dimensional Dirac op- 
erator in quantum field theory. This time, one is interested in extracting the parity- 
violating piece of the effective action Tgff, 

exp (-Teff) = / (PY)(PY) (- I d^x Y(x)/VD,Y(x)) , (3.10) 
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which is obtained by integrating out the fermions. Here, the Dirac operator is mini- 
mally coupled to external fields, such as gauge and/ or metric fields, and, therefore, 
the fermion effective action depends explicitly on them. There is a natural relation 
between the parity anomaly in odd dimensions and the gauge and/ or gravitational 
contribution to axial anomalies in even dimensions [12J. In our case, the imaginary 
part of the effective action that provides the parity violating piece is 

ImFeff = ^J/D / (3.11) 

which in turn is related to the parity violating Chem-Simons action. The derivation is 
based on the formal relation (choosing the branch log(— 1) = in) 

1, det(— A) in ,^^„x 

and it can be made rigorous by zeta-function regularization of the determinant of the 
Dirac operator flSJ. Thus, the Chem-Simons term of topologically massive gauge 
and/ or gravitational theories in three dimensions is induced radiatively |fT9l |20i |2T| 
(but see also [22J for earlier work on the subject). A closely related subject is the struc- 
ture of the induced anomalous vacuum current in odd dimensions and the //-function 
regularization of the vacuum charge of a fermion field (see, for instance, [|23| for a 
general and rigorous discussion of this subject). 



4 Universality of //-invariant in Lif shitz models 

The index of the Dirac-Lifshitz operator i'y^T>^ defined by (|2.1[) (we will use T> in short) 
follows from the integrated form of the axial anomaly. Thus, on a space-time M4 with 
boundaries we have 

lnd{V) = / MR AR) + -1-2 / Me AR)- l^dM,) , (4.1) 

19271^ JMi 19271^ JdMi 2 

in exact analogy with the index of the Dirac operator (|3.4|) , provided that the coeffi- 



cient of the anomalous term in (|2.9|) is the same as for the Dirac operator in a metric 



background. Here, the //-invariant refers to the operator of order z 

zYD; = hy'[D,{-D,,D'' + M^y + {-DkD'' + AfyOi] (4.2) 

acting on two-component spinors on the three-dimensional boundary dM^. 

As before, the Atiyah-Patodi-Singer theorem (|4.1|) implies the following identity 
between //o and the gravitational Chern-Simons action on a compact three-manifold 
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E3 without boundaries 



W, 



csl 



constant . 



(4.3) 



Then, comparison with equation (|3.6|) for the Dirac operator implies that ^70(^3) ought 
to be equal to //d(E3) up to a constant. For E3 ~ S^, in particular, this constant has 
to be zero because S"^ admits a constant curvature metric for which there can be no 
spectral asymmetry for D nor D. Conversely, proving that the //-invariant is universal, 
i.e., 

MS^) = VuiS') (4.4) 

for any given metric on S'^, suffices to show that the coefficient of the anomaly in 
the axial current conservation law (|2.9|) is indeed —l/192n^. This is precisely what 



we will show next by restricting attention to a particular class of metrics that depend 
continuously on a particular modulus (anisotropy parameter). Hence, we will provide 
an alternative derivation of the anomaly by spectral methods, as advertised before. 

The class of metrics that will be used in the following allow for exact computa- 
tion of the jy-invariant. Then, as consequence of the Atiyah-Patodi-Singer theorem, 
equation (|4.4[) should hold true for all metrics on S'^. This way we will also be able 
to provide model geometries for which the index of the four-dimensional fermion op- 
erator can be computed explicitly by spectral flow. Finally, we will derive geometric 
conditions for having non-zero index on certain backgrounds. The presentation we 
follow here has many common elements to our earlier work [7J, but it is now applica- 
ble to Lifshitz operators of arbitrary order z. 

First, to set up the stage, we consider some general aspects of the eigen-value prob- 
lem of the Dirac operator 27' on a three-dimensional manifold E3 with a Riemannian 
metric g acting on two-component spinors. 



zVD,-Y(x) =CY(x); Y(x) 



(4.5) 



where 



(4.6) 



Here, Ei are the components of the inverse dreibeins associated to gij with tangent 
space indices I and co^^i are the components of the corresponding spin connection. 
Also, the gamma matrices 7^ are chosen in terms of the Pauli matrices, as usual. 
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J 






-1/ 



(4.7) 
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and satisfy the anti-commutation relations [7^,7^]+ = 2S^K A useful relationship is 
provided by Lichnerowicz's formula for the square of the Dirac operator |[24|1, 

(zVD,)2 = -D2 + 1r, (4.8) 

which is written in terms of the Bochner Laplacian operator = D/D' acting on 
spinors (rather than scalars) and the Ricci scalar curvature R associated to the three- 
dimensional metric g. 

The eigen-value problem for the more general class of operators 27^ of order 
z = 2a + l, 

zYD/Y(x) = Z Y(x) , (4.9) 

is definitely a more complex problem. Even if the spectrum ^ of the Dirac operator 
can be determined for a given metric g, the eigen-values Z will not be easy to obtain 



in general. Note at this point, based on Lichnerowicz's formula (|4.8|) , that we may 
express the operators Z7'E), in the form 

iyiDi = l^V [d, (^{ij'D^f - + + (^(ij'Dkf -^R + D,] (4.10) 

from which it follows that 27' and 27^ do not commute in general. Their com- 
mutator differs from zero by derivatives of the curvature, which provide a geometric 
obstruction to find a common system of eigen-spinors Y. However, the situation sim- 
plifies drastically when the geometry is homogeneous, since R is the same at all points 
of space and all its derivatives vanish by definition. Only in this case the two oper- 
ators exhibit a common system of eigen-spinors and the corresponding eigen-values 
are simply related by 

Z = ^(^^-Ir + mA . (4.11) 



From now on we restrict attention to homogeneous geometries on Z3 and to be 
more specific we consider the class of Bianchi IX metrics on E3 ~ of the following 
type 



ds^ = j^{a'y+((ry + S^{ay\ , (4.12) 
using the left-invariant one-forms of the group SU{2) that satisfy the relations 

da^ + ^e^jK(r^ A (7^ = . (4.13) 

Here, we are actually assuming that the geometry has enlarged isometry SU{2) xU{l) 
by taking two of the metric coefficients equal to each other. Then, the parameters of 
the metric are the conformal factor 7 and the anisotropy parameter S, which can range 
from to 00. These metrics are often referred in the literature as Berger spheres. 



A\2 , /_2\2 , r2/_3\2 
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The spectrum of the Dirac operator can be completely determined on such geome- 
tries and the same thing also applies to the spectrum of the more general operators 
Z7^D/. If we were considering the more general class of Bianchi IX metrics, without 
imposing axial symmetry, the eigen-value problem of the Dirac operator would have 
not been tractable. For later use we also write the Ricci scalar curvature for this par- 



ticular class of three-dimensional metrics (|4.12|l . 



R = ^{^-S^). (4.14) 

There is a dual basis of vector fields fi associated to cr^ with < fi,cr^ >= ^1 that 
satisfy the commutation relations [//, fj] = —Cij^/k and represent the SU{2) Killing 
vector fields of this particular class of geometries. Both and // can be written in 
terms of three Euler angles, but the explicit expressions will not be needed here. Then, 
the Dirac operator for the class of axially symmetric Bianchi IX geometries (|4.12|l takes 
the following form (see also the Appendix) 



( ifs/S ifi-f2 \ 



\ ifi+fi -ifs/^ J 



+ ^(^' + 2), (4.15) 



where we have taken 7 = 1 without loss of generality. The conformal factor 7 can be 
easily reinstated by noting that 27' Dj scales uniformly as 1/ 1^7, but this will not affect 
the value of the //-invariant on Berger spheres, since it only depends on S. For later 
use we note that the same thing applies to the operators iyiDi that scale uniformly by 
1/ (y^)^'*'^^ and their //-invariant is also be inert to 7 (actually, in this case, we also 
have to rescale appropriately the mass parameter M, but it does not really matter since 
the //-invariant will turn out to be independent of M as well). 

The spectrum of the Dirac operator (|4.15|) on Berger spheres is given by the follow- 
ing discrete set of eigen-values [|25il (but see also [|261 ) 



C± = I ± ^ ^^S^pci + ip-c,)^ ; (p, c,) G N2 (4.16) 

with multiplicities p + cj for each pair {p,ci)- The positive integers p and q are not 
ordered, which means that the conjugate pairs {q, p) also yield the same eigen-values 
^± with the same multiplicities. There are additional eigen-values arising formally 
for q = 0, 

but their multiplicity is 2p rather than p. Thus, the complete spectrum of the Dirac 
operator is ^± and ^0 forming a two-dimensional state lattice which is labeled by the 
quantum numbers {p,q)- The eigen-values ^+ and ^0 positive definite for all 5, 



12 



whereas ^_ can assume positive as well as negative values by tuning their parameters. 
Likewise, the spectrum of the more general operator 27' takes the form 



Z = c(^^^ + l{S^-4:)+M^j (4.18) 

following from (|4.11|) and (|4.14|) with 7 = 1. The corresponding eigen-values Z± and 
Zq are obtained by setting ^ = ^± and ^0 respectively. 

We also note that the sign of the eigenvalues ^ and Z is the same irrespective of a. 
and M, since 

^^ + hs^-4:)>0, V^G[0,oo). (4.19) 
8 

This can be verified directly for ^ = ^± or ^0 arid it is important for computing the 
^/-invariant of the corresponding operators by separating the contribution of their pos- 
itive and negative modes. Although ^+ and ^0 are positive definite for all values of 
S, I,- can be either negative or positive depending upon 3. The eigen-values ^- can 
also become zero at special values of 5 for which there exist positive integer solutions 
{j),c\) to the equation 

8^ = 2^ ApqS^ + {p - . (4.20) 

The corresponding zero modes Y are often referred as harmonic spinors in the math- 
ematics literature ||24ll25ll , and, clearly, they can only exist for S > A. The presence of 
zero modes should be properly accounted in the computation of the //-invariant. 

There are some important points that need to be emphasized before proceeding 
further. The anisotropy parameter 3 accounts for the distortion of the three-sphere. 
The value ^ = 1 corresponds to the round metric on S^. When S decreases the sphere 
is squashed until it becomes fully flattened in the extreme case ^ = 0; in the latter 
case collapses to but the geometry is completely regular there. When S increases 
the sphere is stretched until it becomes singular as — )■ 00. Note that the Ricci scalar 
curvature (|4.14|) is positive definite for S < 2 and it turns negative for S > 2. Then, 



by Lichnerowicz's formula (|4.8|) , the Dirac operator cannot exhibit any zero modes for 
S < 2. Zero modes are allowed to exist only for S > 2. Note at this point that when 
the curvature of space is positive definite, all eigen-values ^_ are negative and clearly 
the same thing applies to all eigen-values Z_ . For negative curvature, however, level 
crossing becomes possible provided that ^ > 4. The critical value (5 = 4 is the thresh- 



old for the existence of solutions to equation (|4.20|) . This value is far beyond Lich- 
nerowicz's bound S > 2, and, thus, in this context, it corresponds to Berger spheres 
with sufficiently negative curvature. At 5 = 4: the eigen-value ^_ with p = q = 1 
becomes zero and then it turns positive for all ^ > 4. As 3 increases further, there 
are more and more eigen-values ^- that become positive, until they all undergo level 
crossing in the extreme limit S ^ 00. 

The computation of the //-invariant fjniS^) relies on the separation of positive and 
negative modes. First, we assume that ^ < 4 so that the eigen-values Z_ are all neg- 
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ative, whereas Z_|- and Zq are all positive. Later we will extend the result to Berger 
spheres with ^ > 4 by counting the number of modes that have undergone level 
crossing. Thus, according to definition. 



//D = lim7D(s) = lim V (signZ)|Z| 

eigenvalues 

we have to evaluate the following sums when S < 4, 

'5^ 



(4.21) 



p,q>0 



p>0 



+ p 



+ X 



5^ 



2 x2 



V + X) + V(^^-4)+4^2^2 



" 2 



+ P) +^(^^-4)+4^2^2 



+ 



(4.22) 



setting X = ^y^S'^pq + (p — for notational convenience. The first line refers to 
the contribution of Z+, the second to Z_ and the third to Zq. Note that in writing 
down (|4.22|) all eigen-values ^ in (|4.18|) have been rescaled by a factor of 23, which is 
irrelevant, however, for the final result that is obtained at s = 0. 

Next, we compute the individual sums by expanding all fractions in powers of 5 
(up to the appropriate order) and then set s = 0. The contribution of the modes Zq{S) 
is the easiest to evaluate. Using the power series expansion 



p(2a+l)s 
+ • • 



,2/^2a + l 2a(2M2-l)\ sS^ f , ,,0 ^.n 

1 - 3(^2 I — + ^ ^ L ) + _ ( (s + i)(2a + 1)2 - 2a(2a + 3) 



2p 



8p2 



(4.23) 



we find that the last group of terms in (|4.22|) assumes the following expansion written 
in terms of Riemann zeta-functions. 



2^(2^ + 1)5-1 



{2a + l)C{{2oi + l)s) + 4a(2M2 - l)C((2a + l)s + 1^ 



+ - 



s + 1)(2^^ + ly - 2(x{2(x + 3) ^((2a + l)s + 1) + 



(4.24) 



The terms that are omitted vanish at s = as they contain the factor s^{{2oi + l)s + n) 
with integer n > 2 and ^(s) is absolutely convergent for Re s > 1; they include a 
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term of the form S^/p'^ in (|4.23|) as well as all terms of order and higher which come 
multiplied with 1/ with n > 2. Then, taking into account that — 1) = —1/12, 
^(0) = —1/2 and that ^(s) has a simple pole at s = 1 with residue 1 (and, thus, 
(2a + 1)s^((2a: + l)s + 1) equals 1 at s = 0), we obtain 



4a 



2a: -1 



(4.25) 



The contribution of the modes Z±{S) to the //-invariant is more difficult to extract 
because of the double sums that are involved. Again, we expand the fractions in 
power series of S 



5^ 



±— + X) + — (^2 _4) ^4^2^2 



2 



X(2a+l)s 



2X 



X2 



^^4 |^ (s + l)(2a + 1)2 -2^(2^ + 3) ^ a[(s + l)(2a + l) - (2a-l)](2M2-l)) - ^ ^ 



8X2 



X3 



^3 {2cc{2cc - l)(2a + 7) + (s + l)(2a + l)[(s + 2)(2^^ + 1)2 - 6a (2a + 3)] 



+ 



(4.26) 



omitting all terms of the form l/X^+2 with n > 2 that do not contribute to the final 
result when s = 0, as will be seen shortly. Then, the first two sums in the general 



expression (|4.22|l for >/d(s), which we denote respectively by l±{s), are combined to- 
gether as follows. 



7+(s)-7_(s) 



2s/ 

48 
+ • 



^ |a[2a - 1 - (s + l)(2a + 1)](2M2 - l),^* + 



(2a + l)s + 3 



(2a(2a - l)(2a + 7) + (s + l)(2a + l)[(s + 2) (2a + 1)2 - 6a(2a + 3)]) } 

(4.27) 



setting for convenience 



m = E 



p + q 



P + q 



21s 



(4.28) 
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The function /(s) is absolutely convergent on the complex s-plane for Re s > 3/2, 
which justifies the suppression of the higher order terms in the power series expansion 
above (all such terms are multiplied with s and vanish at s = 0). Furthermore, the 
residues of the function /(s) at the two special points s = 1/2 and 3/2, which are 
relevant for the calculation, are (see, for instance, ||27|) 

5^-1 1 
res/(s)|s=i/2 = , res/(s)|s=3/2 = ^ , (4.29) 

and, therefore, the resulting contribution of these terms to the //-invariant is 

( ^ Vloc 2 (- 10a -1 



/,{0) - /_(0) = - ^1 + — (2M^ - 1) j - 3 - 4{2^Tr I • <*-3°) 

Putting all together, it turns out that the //-invariant is independent of the exponent 
OL and the scale parameter M. The final result reads 

7d(0) = J+ - J_ + Jo = -^(1 - d'f. (4.31) 

As such, it coincides with the //-invariant of the Dirac operator on Berger spheres with 
anisotropy parameter 5 and it is related to the gravitational Chem-Simons action as 
follows (see also the Appendix), 

rWcs (Berger) = 1(2 - 2^^ ^ 5^) = -^//d(0) + ^ , (4.32) 



in agreement with formula (|3.6|) for the Dirac operator. This proves our claim (|4.4|) for 



the //-invariant on Berger spheres with ^ < 4 and suffices to show that the coefficient 
of the anomaly is indeed universal. 

Next, we examine the case ^ > 4 that involves level crossing. Every time an eigen- 
value Z_ crosses from negative to positive values, //o jumps by 2 to //id + 2, since the 
spectral asymmetry, as calculated for S < 4, changes by 2. Let us denote by 



C{Sc) = {ip,q) e N^; 3i = 2sJ ^VH^l + (p - qf] (4.33) 

the set of positive integers {p,q) that account for harmonic sptnors at each one of 
the special values < ^ arising for fixed ^ > 4. Since the zero modes at 5c have 
multiplicity p + q, the quantity 

S(^)=Ef E (P + ^)l (4-34) 

Sc<S (p,q)eC(Sc) 

provides the total number of modes that crossed from negative to positive values as S 
was varying from ^ < 4 to any given value ^ > 4. Thus, for ^ > 4, the //-invariant of 
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the operator D on Berger spheres is shifted by twice the number of these modes and 
equals 



S{S) is also inert to the parameters a. and M which shows that level crossing is the same 
as for the Dirac operator. It can be easily seen that S{S) is always an even integer, but 
there is no closed expression for it for general values of S. 

The results we have presented here generalize our earlier work on the subject fZl to 
all values of the anisotropy scaling parameter z and show that the axial anomaly is an 
infra-red phenomenon, as in relativistic quantum field theory The universality of the 
//-invariant also shows that the parity breaking piece of the induced fermion action 
in three-dimensions is provided by the Chem-Simons term with the same coefficient 
as in (|3.11|) irrespective of the order of the fermion operator Z7^D,. Finally, the model 
geometries that were used to perform the calculation provide simple examples for 
which the index of the four-dimensional fermion operator can be computed explicitly 
by spectral flow as Ind(P) = ±S(^) (the sign depends on the choice of orientation). 
In the next section, we consider four-dimensional backgrounds that can give rise to 
violations of the chiral charge conservation law in gravitational theories, and, thus, 
lead to baryon and lepton number violation when the index differs from zero (note at 
this end that the axial anomaly will not obstruct the conservation law AQ5 = if the 
index is zero). 

5 Applications to non-relativistic theories of gravity 

So far the geometry has been treated as background for fermion propagation without 
reference to any metric field equations. Space-times of the form I x E3 with Euclidean 
line element (|3.7|) will lead to violations of the chiral charge conservation if the gravi- 
tational field equations for the metric gij{t, x) induce level crossing for the fermion op- 
erator on E3 as t varies in I. In Einstein gravity this cannot happen for various reasons 
(see, for instance, ||28|), but in non-relativistic gravitational theories of Hof ava-Lifshitz 
type ||29| there are regular solutions that can support such exotic effects. 

For definiteness, we concentrate on gravitational instanton solutions in 3 + 1 di- 
mensions with SU{2) X U{1) isometry for which E3 ~ S^^ is a Berger sphere with 
metric (|4.12|) . In Einstein gravity these are the familiar Taub-NUT and Eguchi-Hanson 
instantons for which I is the semi-infinite line and the anisotropy parameter S changes 
between and 1 as i varies between the two end-points of I; in these cases one side of 
I terminates at a removable singularity, which is a nut or a bolt respectively. The field 
equations do not allow S to exceed 1 without ruining the regularity and completeness 
of the four-dimensional metrics. Thus, in this case, level crossing cannot occur and 
the index of the four-dimensional Dirac operator on such metric backgrounds is zero. 
More generally, for non-compact four-metrics with non-negative Ricci scalar curva- 




(4.35) 
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ture, there can be no bound state solutions of the Dirac equation; if such states existed, 
they would have been covariantly constant, and, hence, non-normalizable | |24| leading 
to contradiction. 

On the other hand, in Horava-Lifshitz theory with detailed balance (meaning that 
the potential term of the action is derivable from a local superpotential functional 
W[g]) the instantons are eternal solution of the gradient flow equation for the three- 
metric g derived from W ||30| . In this case, J ~ K and the instanton interpolates 
between two distinct vacua of the three-dimensional action W (viewed as fixed points 
of the gradient flow equation) as t ranges from — oo to +00. Typically, on S^, W[g] 
admits a maximally symmetric vacuum for which S = 1, but there should be other 
less symmetric vacua that coexist and provide the other end-point of the instantons. 
Then, within our Bianchi IX mini super-space model, we pose the following problem: 
under what conditions W[g] admits another less symmetric vacuum with ^ > 4, in 
order to have violation of chiral charge conservation induced by instantons, and what 
is the phase of the non-relativistic gravitational theory that admits such instantons as 
bona fide Euclidean solutions? The results we describe below show that Q5 cannot be 
always conserved in Hof ava-Lifshitz gravity with anisotropy scaling exponent z > 3. 
For lower values of z there seems to be no window in the space of couplings for having 
non-conservation of chiral charge. Curiously, Hof ava-Lifshitz gravity becomes power 
counting renormalizable only for z > 3, which was the reason to view it as substitute 
of general relativity in the deep ultra-violet regime in the first place [29|. Thus, it is 
interesting to study the problem of chiral symmetry breaking in these models. 

Recall that the superpotential functional W[g] is taken to be the action of three- 
dimensional gravity on Z3 augmented with higher derivative terms. A simple choice 
is provided by topologically massive gravity [22] containing the gravitational Chern- 
Simons term (|3.5|) and a three-dimensional cosmological constant A^, 



VMG[g] = 4- / d'x,/d^ {R - 2Aw) + ^Wcslg] . (5.1) 

The corresponding theory of Hof ava-Lifshitz gravity has anisotropy scaling parame- 
ter z = 3 unless the Chern-Simons term is absent in which case z = 2. The vacua of 
topologically massive gravity satisfy the following field equations, which are split on 
purpose into traceless and trace parts. 



^-\Hti-—Cii = ^ and J^ = 6Aw, (5.2) 

where C,y is the Cotton tensor of the three-metric g. 

For the class of Bianchi IX metrics (|4.12|) , the traceless part of the field equations 



reduce to a single algebraic relation among the anisotropy parameter 5 and the con- 
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formal factor 7, 

Thus, apart from the maximally symmetric solution 3 = 1, there is one more solution 
for 

l + AL = (5.4) 

provided that the Chern-Simons coupling 00 is negative (for a given choice of orien- 
tation on S^). The two metrics can coexist as vacua of topologically massive gravity 
only when is non-negative. This follows from the trace part of the field equations 
K = 6Aw using (|4.14|l . Then, the corresponding instanton of Horava-Lifshitz grav- 



ity interpolates between two Berger spheres with the same non-negative curvature, 
and, thus, there can be no net level crossing in the spectrum of the three-dimensional 
fermion operator. The only way to have level crossing is to be in a unimodular phase 
of the theory in which the trace part of the field equations decouples from the dynam- 
ics and the two end-points are allowed to have different curvature (positive for ^ = 1 
and sufficiently negative for the other solution (|5.4|) ). 

This possibility is naturally realized when the generalized DeWitt metric in super- 
space arising in the canonical formulation of the theory [29J 

has A = 1/3 and projects any tensor to its traceless part (Gjjki has similar form and 
properties). Then, the volume of space is preserved by the dynamics in 3 + 1 dimen- 
sions. Since the volume of the Berger sphere (|4.12[) is fixed Vo^S"^) = IGti'^Sj^^^, we 
can rewrite the defining equation of the anisotropic fixed point (|5.4|) as follows by 
eliminating 7, 

(2 \ 3 
-^J . (5.6) 

Thus, in this case, the index of the Dirac-Lifshitz operator will be non-zero on the 
corresponding gravitational instanton background provided that ^ > 4, which in turn 
imposes a lower bound on the volume of space written in terms of the ratio of the 
couplings to CO as follows. 



Vol(S^) > 138247r2 ( -— ) , (5.7) 



CO 



reproducing the results reported in f7\. Note that the second fixed point (|5.4[) is absent 
when there is no Chem-Simons term in and, hence, there is no instanton to 

induce violation of chiral charge conservation when the anisotropy scaling exponent 
is lowered to z = 2. 
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Next, we consider the case of z = 4 Hof ava-Lifshitz gravity in 3 + 1 dimensions 
(see, for instance, 113 TP choosing as W[g] the action of new massive gravity on E3, [32], 



WNMG[g] = 4- / d^xy^ {R - 2Aw) + ^ / dhy^ (R.jR'j - IrA . (5.8) 



We want to examine if the results described above persist in this case, leaving again a 
window for violation of chiral charge conservation by gravitational instanton effects. 
The choice (|5.8|) is not the most general one that can be made. One may augment it 
with a gravitational Chem-Simons term, as will be discussed briefly later, and also 
change the relative coefficient of the quadratic curvature terms in (|5.8|) , which, how- 
ever, will not be discussed here for simplicity (such a generalization spoils some nice 
properties of new massive gravity as higher curvature extension of three dimensional 
Einstein gravity, but it is legitimate in the context of z = 4 Horava-Lifshitz gravity). 

For Berger spheres, the traceless part of the classical equations of motion of new 



massive gravity reduce consistently to the following algebraic relation (see [|33| for 
details). 

Thus, apart from the maximally symmetric solution ^ = 1 we also have another solu- 
tion for 

i + ?%:|i?!)=o. (5.10) 

Since we are ultimately interested in the coexistence of vacua with 5 > ^, we impose 
the restriction > on the coupling of the quadratic curvature terms (otherwise, if 
5^ < 4/21, the existence of a second solution will require < 0). The trace part of 
the field equations yields 

^ = -^Aw + ^(3<S^ - 4){7^^ - 4) . (5.11) 



Setting S = 1, we find that Aw is non-negative for > 0. The other solution (|5.10|) 



does not satisfy the trace relation for ^ > 4, since the left-hand side will be negative 
and the right-hand side positive. Thus, the only way that both solutions can coexist 
when (J > 4 is to decouple the trace part of the field equations from the dynamics. 

As before, this only happens in the unimodular phase of the theory corresponding 
to the choice A = 1/3 in the DeWitt metric. Then, the volume of space remains fixed 
and can be expressed in terms of 3 by eliminating 7 from the fixed point (|5.10|) . We 
have, in particular, the following relation 

Vol(S3) = ^^(21^2 _ 4)3/2 > 8^2(83)3/2^ ^ (5 -^2) 
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using ^ > 4 to establish the lower bound in the volume of space in terms of the cou- 
plings. This inequality provides the window for having violation of chiral charge 
conservation in the corresponding 3 + 1 theory of Horava-Lifshitz gravity. 

Further generalizations arise by adding the gravitational Chern-Simons term to 
Wnmg is] yield the so called three-dimensional generalized new massive gravity 
[32J. This theory has a maximally symmetric solution with ^ = 1. Other Berger sphere 
solutions satisfy the traceless part of the field equations, which boil down to 

1 k^(4-21^ 

3 + 2^+ 48m^j =°- ^^-^^^ 

For simplicity, we assume > so that there is only one admissible solution of 
equation (|5.13|) , whereas co can take either positive or negative values. The trace part of 
the field equations is the same as before, (|5.11|) , because the Cotton tensor is traceless. 
Thus, for > 0, (|5.11|) cannot be possibly satisfied for both ^ = 1 and (5 > 4 and we 
are forced again to consider the unimodular variant of the theory that corresponds to 
A = 1/3. Then, the volume of space is preserved throughout the evolution and it is 
given by 



Vol(S' 




> 6471^ + w + rp^ . (5.14) 



The lower bound is obtained by demanding (5 > 4 so that chiral charge conservation 
can be violated in the corresponding Horava-Lifshitz gravity and it encompasses the 
two simpler cases that were discussed before. The analysis is more intricate when 
< 0, as there can be more than one vacua with S ^ 1, and, hence, more instantons 
that interpolate between them; the details are left to the interested reader. 

Similar considerations apply to theories based on Born-Infeld generalizations of 
three-dimensional gravity [34] that include higher and higher derivative terms. The 
anisotropy scaling exponent z of the corresponding Horava-Lifshitz theory increases 
accordingly and there can be several Berger sphere solutions that coexist and support 
the end-points of instantons. We expect that the general features of our results will re- 
main qualitatively the same in these cases too, leaving windows in the space of grav- 
itational coupling parameters for having chiral symmetry breaking by gravitational 
instanton effects. This novelty differentiates the non-relativistic gravitational theories 
of Lifshitz type from Einstein gravity and will have phenomenological consequences 
in bouncing models of the early universe if Horava-Lifshitz gravity plays a role. We 
also expect these results to be valid more generally beyond the mini-superspace model 



21 



geometries, which were only used here for illustrative purposes. 

In all examples that were considered above, it was found that non-conservation of 
chiral charge can only occur in the phase of Horava-Lif shitz gravity with A = 1/3, 
where the volume of space decouples from the dynamics. This phase of the theory is 
not inflicted with many of the pathologies that otherwise haunt non-relativistic mod- 
els of gravitation. It may very well be that only this phase becomes relevant in the 
deep ultra-violet regime, as substitute for ordinary gravity, and not its variants with 
A ^ 1/3. Flowing from the ultra-violet to the infra-red is a difficult problem that could 
be properly investigated within the asymptotic safety scenario for quantum gravity 
(for a review see, for instance, [35] and references therein) treating A as a coupling that 
undergoes renormalization. Proper treatment of the problem may also require turning 
on other operators that are lying outside the class of Horava-Lifshitz models to ensure 
well-behaved transition to ordinary large scale physics. All these issues are currently 
under consideration. 



6 Conclusions 

We have investigated the occurrence of axial anomalies in non-relativistic fermion 
theories of Lif shitz type which are minimally coupled to gauge and /or gravitational 
fields in 3 + 1 space-time dimensions. Our main conclusion is that the results are iden- 
tical to the relativistic case, hereby generalizing earlier work on subject HI |Z1 to theo- 
ries with arbitrary anisotropy scaling exponents. This is not surprising in retrospect 
conforming with the general idea that the axial anomaly (unlike others) is an infra-red 
phenomenon. Also, the universal form of the ^/-invariant can be intuitively under- 
stood (prior to zeta-function regularization) by noting that the operator (— + M^)* 
is positive definite; hence, the number of positive and negative modes is the same for 
both operators 27' D, and 27^ D,. It was confirmed that the same conclusion holds for 
the regulated 7/ -invariants. 

Although the result has its own value, we have also considered some applications 
in the context of non-relativistic gravitational theories of Lifshitz type and found that 
(unlike ordinary gravity) violation of chiral charge conservation can be induced by 
instantons - and possibly other configurations - in certain regions of their parameter 
space. 

There are a few technical questions that remain open and should be addressed 
carefully elsewhere. It is not yet clear whether there are any restrictions on the al- 
lowed values of the anisotropy scaling parameter z that insure self-adjointness of the 
Dirac-Lifshitz operator in the strict mathematical sense. Also, it is interesting to gen- 
eralize our results to higher dimensions and compare with the anomalies of the corre- 
sponding relativistic quantum field theories. Finally, it remains to explore the role of 
anomalies in condensed matter physics systems that exhibit quantum criticality. 
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A Chern-Simons action for Berger spheres 

In this appendix we provide the components of the connection one-form of the Berger 
sphere with metric ds^ = 7 ((c7^)^ + (c^^)^ + S^{a^)^). These expressions are needed 
for writing down the Dirac operator on such backgrounds as well as for evaluating 
the gravitational Chem-Simons action Wcs (Berger) used in the main text. We have, in 
particular, 

w\ = ^-^(^ , 0)^3 = , ^^^3 = (A.l) 

noting that they are all independent of the scale factor 7. 

The gravitational Chern-Simons action can be easily evaluated in steps as follows, 

Wcs = Tr ^o; A i/o; + ^o; A a; A o;^ 

= 2 J^^ (^co^2 A dco^i + io^3 A dco^i + A dco^2 + 2co^2 A (0^3 A co^i^ 

= ((5^-2(52+2) / a^Aa^Aa^, (A.2) 

using along the way dcr^ + cr^ Acr^ = and cyclic permutations. Since the integral 
over of the volume form Aa^ A cr^ equals 1671^, we obtain the final result 

Wcs (Berger) = 16n^ {S^ - 2S^ + 2) . (A.3) 

The extrema of this action, as function of ^ G [0, 00), occur at ^ = 1 and ^ = in 
accordance with the fact that the round sphere and the fully squashed Berger sphere 
are the only metrics in this class that are conformally flat. The minimum corresponds 
toS — 1, whereas (5 = is a local maximum. 
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